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We are interested in solving the coupled set of equations 




( 1 ) 



For the problems we are interested in, these equations do not reduce to any 
known special functions. These equations arise in the study of combat models 
where x and y are opposing forces and a(t) and b(t) are the attrition 
rate coefficients related to the effectiveness of weapons. Lanchester [1] 
first looked at these types of models for aircraft battles, using constant 
coefficients (in which case the solution is in terms of exponentials) and so 
these types of equations are referred to as Lanchester equations. There are 
also special cases of the coefficients for which the solutions to (1) are 
expressible in terms of the generalized Airy functions [2]. 

We are interested in problems where the coefficients a(t) and b(t) 
are positive and continuous. In that case we expect that we can find solutions 
which are linear combinations of monotonic, exponential like functions. If 
we have such a set, then the analysis of the forces necessary to guarantee 
a win for one side, say x , is very much simplified. So we are interested 
in determining exponential like solutions for (1), and error bounds for the 
solutions. 

We convert (1) into a single second order differential equations for x(t) 




( 2 ) 



x(0) = Xg , a-’(O) ^ (0) = -y^ 



( 3 ) 



Equation (2) is still awkward. Let us define two useful quantities, the 
intensity of combat I(t) and relative fire effectiveness R(t) , 



I(t) ^ /a(t)b(t) 



( 3 ) 



R(t) H 




(4) 



A change of independent variables will simplify (2). Let 



then (2) becomes 



- / 

0 



I(s)ds 



9 





dR dx 

“ 



X = 0 . 



(5) 

( 6 ) 



We now reduce (6) to the standard Liouvil 1 e-Green form by the usual change 
of dependent variable 

H - . (7) 

/r 



So (6) becomes 



d^W 

dx2 



1 + /R 



dx 



(1//R) 



W = 0 



( 8 ) 



Then the standard estimates for (8) say that there are two solutions to (8) 



W = c-je"^ + 026^ , 

i.e. there are two solutions to (2) 

t t 

( / I(s)ds -/ I(s)ds I 

What is the error in using (9)? A number of error estimates are avail abl 
[^J» [4]- It turns out that they are not quite good enough for our purposes. 
The estimates depend upon whether one is interested in the dominant solution 
or in the subdominant solution. 



( 4 ) 



Looking first at the dominant solution, , we take 



= e”^[l + h^(x)J 



( 10 ) 



Then the error h-jCi) satisfies 



wnere 



I t I 

h^ + 2h-| = il- (1 + h^) 

.2 

,p = /r-2- ( 1 //^ . 

dx 



(11) 



( 12 ) 



Converting (11) to an integral equation with zero initial conditions we have 

H,(t) = i/'[l - t h,(s)Jds . (13) 

Our error estimates are based on Wil let's generalization of Gronwall's lemma. 
Lemma 1 . Let 



and let 



then 



h(x) = / K(s,x)>J'(s)[J(s) + h(s)]ds 
0 



|K(s,x)| < Q(s)P(x) , 



|h(x)| < P(x) / I J(s) |Q(s) l’l'(s) I exp|/ P(z)Q(z) |i(i(z) |dz|ds . (14) 

Proof: See Wil let [5]. 

An immediate corollary is: 

Corol 1 ary . Let M(z) = Q(z) max { | J(z) | ,P(z)} then 



|h(x)| < P(x) 



expj/ M(s) |ij;(s) |ds| - 1 



Proof: Integrate (14), using M(z) . 

In the case of equation (13) 

K(s,x) ~^(1 “ 6 

J(s) = 1 . 



(15) 

OED 



- 1 M ,-2(x - s)j ^ 1 _ g-2xj 



( 5 ) 



-2t 1 

Thus we can take P(x) = 1 - e and Q(s) = j . Then lemma 1 gives 
Theorem 1 . The dominant solution to (8) satisfies 



W-| = e^[l + h-| (t)J 



< e 



1 + (1 - e"^^)(exp|/ ^ U(s)|ds| - 1) 



As for the subdominant solution, we take 



and then satisfies 



which converts to 



W 2 = e" [1 + h 2 (x)] 

I I I 






" 1 r^ZCx-s) 



- 1 



4 »(s)[l + h 2 (s)]ds 



Now we have 



-2s 



< e 



-T 



(16) 

QED 



(17) 



(18) 



K(s,x) = 1 e2(^-s). 1 < ]- (e2^ - l)e 

Then we have 

Theorem 2. The subdominant solution to (8) satisfies 



W 2 = e"^[l + h2(x)] 

1 + (e^"^ - 1) / / (1 " 

0 ^ s ^ . 

Observe that we did not use the corollary on (19). 

We can see that our estimates depend very strongly on ip , given by (12). 

For the Lanchester equations (1), the coefficients a(t) and b(t) , for 

many applications, can be expressed by 



. (19) 

QED 



a(t) = k,(t + c)^ 

a 

b(t) = kj^(t + c + a)' 



(20) 



( 6 ) 



where c is the "starting" parameter, and A is an "off-set" parameter 
(so y can have a different firing range than x ). These are referred to 
as the "power rate" coefficients. For future use, let 



B = ATiT 
a b 



r _ V + u + 2 
5 ^ 



We first consider the case of no offset, i.e. A = 0 . 

t 

T = B / (t + c) ^ dt 
0 






and 






u - u (u + 3yi + 4) 
i 6 / . , „6 \ 2 

(6t + C ) 



We see that for -1 < u < u that ^ is negative, and 

- 1 



exp / I ds 
0 



= -1 + exp 



u - p 



326 B 






Then (see (5)) 



exp - 



Let 



'u - pVu + 3p + 4 



+ 6 X 



325B A,6 



= ^ ' 326B^^ (u + 3p + 4) 



and = exp {y/^^) 



Then (16) becomes 



w^ < e^( exp (~y/(A‘^ + 6x)) + e"^^ - D^e"^^ exp {-y/{A^ + 6x)) 



-2x 



( 21 ) 



( 22 ) 



(23) 



as an error term for the dominant solution. 

For the subdominant solution the corollary gives too crude an estimate. 
Even with using the lemma, the error estimate is not good. Integrating by 
parts we get 



( 7 ) 



= e“ [1 + t)J 

< e" ^ [1 + 1)^1 - exp|^/ (1 - e'^^)||-|dzj 

+ / lf|e^p[/ 0 - 
0 ^ Ls 



dzlds). 



Working this out we get 



W 2 ^ e 






(25) 



This is not a very exciting bound, but we are unable to do any better, using 
the Willet result. 

We can, however, alter the Willet result. The subdominant solution error 
also satisfies the equation 

hj = j/*[l - e-2(= - + h2(s)]ds . (26) 

T 

We now prove 

Lemma 2 . If h( x) satisfies 

h(x) < / K(s, x)i^(s)[J(s) + h(s)]ds (27) 

T 

and 

lK(s,x)| < P(x)Q(s) (28) 

where K, tp, and J are all s 0 , then 



h(x) < P(x) 



/ iJ;(s)Q(s)J(s)exp / P(a)Q(a)<l^(a)da ds 



(29) , 

I 



Proof . Eq. (27) can be written, using (28) 

pW < / Q(s)iJ^(s)J(s)ds + / Q(s)4«(s)h(s)ds . 

' X X 

Differentiating 



( 8 ) 



(^) ^ -Q(t)*(t)J(t) - q(,)*.(t)h(t) ^ 

Therefore 

t 

(^) + QtP(|) ^ . 

Integrating we get 

2 / q(s);(/(s)J(s)expj^/ P(a)Q(a)iJ;(a)da^ds 

+ K exp^/ P(a)Q(a)i(j(a)da^ . 

But h(») = 0 , so K = 0 . 

For equation (26) we have 



K(s,t) = - e~2(s - x)J < 1 . ^ ^ 



QED 



Thus 



hp < / 






r S 



exp 



= exp !/ l|2idc,[ - 1 



/ j (a)da 



ds 



(30) 



This is the same bound found by Olver [3]. 

For the case of the power attrition coefficients with no offset we have 



which is a fairly reasonable bound. 

For the case where p = u = 1 , the linear case, we can also get some 
explicit results with offset. That is 



Then 



a(t) = k^(t + c) 
b(t) = k|^(t + c + A) 

, = .^)j 



(32) 



( 9 ) 



where 



, , 1 . 

Then 

1 ^ 

^ / t(a)da 

■^0 

^ _J U(t + c)^ - 6A(t + c)^ - 12A^(t + c) - 7A^ 

16BA^ 1 C(t + c)(t + c + 



4c^ - 6Ac^ - 12A^c - 7A^ j 
rC(C + A)]^/^ ) 

while 

oo 

J / 4'(a)da 

T 

^ _J L _ 4(t + c)^ - 6A(t + c)^ - 12A^(t + c) - 7A^ 

16BA^ 1 [(t + c)(t + c + A)]^/^ 

Both integrals are monotone increasing. 

For the general case 



(33) 



(34) I 

I 

t 

I 



J / 4-(s)ds 
^0 

1 n ( (p + |p^)(n^ + An + A^/4) - (u + u^/ 4 )(n^ - An + A^/4 

2/R^ (q < (7"^~A/27~^^^ I 



1 

2 



(n - A/2)^ 



(n + A/2) 



1 + u/2 



dn 



where n = t + c + A/2 . 



(35) 



( 10 ) 
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